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Abstract 

The absorbing "ricepile" model with stochastic toppling rules has 
been numerically studied. Local limited, local unlimited, nonlocal lim- 
ited and nonlocal unlimited versions of the absorbing model have been 
investigated. Transport properties and different dynamical regimes of 
all of the models have been analysed, from the point of view of self or- 
ganized criticality (SOC). Phase transitions between different dynam- 
ical regimes were studied in detail. It was shown, that the absorbing 
models belong to two different universality classes. 
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1 Introduction 



Self organized criticality (SOC) is widely studied phenomenon in the last 
ten years. Theory of SOC, proposed by Bak, Tang and Wiesenfeld 0] de- 
scribes dynamical behaviour of many particle systems with local interactions. 
Paradigmaticaly, the description is based on the dynamics of a pile of sand. 
If the sandpile is randomly driven by slow addition of sandgrains, the slope 
of pile grows up and after some time local stability conditions are violated 
somewhere on the pile surface. The avalanche starts to slide down the slope. 
This type of dynamics is easily modeled by a cellular automaton. Such a 
model "sandpile" is defined on a large n- dimensional lattice. The avalanche 
dynamics is, under the action of slow drive, governed by the local critical 
conditions (such as the local critical slope, for example), and local toppling 
rules. This dynamics leads to the steady state called self organized critical 
state, characterized by the critical scaling of the avalanche size distribution: 

p{s,L) = s-^f{^) (1) 

In (|I]), L is the system size and s exhibits the size of the avalanche. Crit- 
ical exponents r and D, depend on the particular model. Different models 
can be divided into the various universality classes, defined by a specific set 
of the critical exponents |10|, ||] . 

A natural step from the sandpile model systems, leads to the investigation 
of real piles of granular material from the point of view of SOC dynamics. 
Several efforts have been made in this direction ||^, ^, but no clear evi- 
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dence of self organized critical dynamics (|T]) has been found. Finaly, in 1996 
an experiment has been done by the group of experimentalists and theoretists 
in Oslo 0]. In the Oslo experiment dynamical behaviour of the driven qua- 
sionedimensional pile of rice has been investigated. The avalanche sizes in 
the steady state were measured in terms of dissipated potential energy. Two 
type of grains were used (elongated and round ones), showing completely 
different dynamics. In the case of the ricepile, consisting of elongated grains, 
SOC state has been established, in which the avalanche size distribution has 
a power law character with critical exponent r ^ 2.02 . 

Soon after the experimental results were published, the model "ricepile" 
cellular automata were suggested and numerically studied 0, §], .The 
model ricepile is, in principle, cellular automaton defined on onedimensional 
lattice, with randomness incorporated into the toppling rules and with de- 
terministic drive. Changes in toppling rules are often manifested by different 
dynamical behaviour of the model and different universality class into which 
the model belongs, |T0[ 

In this paper we study a model, which exhibits a modification of the 
two threshold ricepile model 0]. In ||, 0] the gravity effects, grain fric- 
tion and the local conditions on the pile are described by the parameter p. 
Two thresholds, namely the critical threshold and the gravity threshold are 
defined. Thresholds governs the movement of the grain on the pile surface. 

We removed the second, gravity threshold of the ricepile model. This 
causes, that the model, depending on the parameter p value, has absorbing 
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properties and interesting dynamical regimes [|TT|. One threshold model in 
two dimensions has been studied by Tadic and Dhar |jl2 . 

Here, different versions of the absorbing model are defined using differ- 
ent toppling rules. Dynamical behaviour of the local limited (LLIM), local 
unlimited (LUNLIM), nonlocal limited (NLIM) and nonlocal unlimited 
(NUNLIM) absorbing model is studied and two distinct universality classes 
are recognized. 

2 Ricepile models 

The experimental results of the Oslo group motivated further theoretical 



studies of the avalanches and dynamics of granular material [0, |TJ] . The 
main question is, which physical properties of the pile granular material, are 
important for SOC state to be established. What is, for example, the role 
of friction, what is the role of the grain shape, the gravity and inertia of 
rolling grains? The ricepile experiment reveals, that no SOC is possible if 
the ricepile consists of round ricegrains. On the contrary, the dynamics of 
pile consisting of elongated grains, is self organized critical ||^ . Certainly, 
the shape of grains, and the additional effects related to the shape, such 
as better packing of grains due to the elongated shape, supressed rolling of 



grains and thus supressed inertia effects, are of great importance p|, [14 



Ricepile models are cellular automata, in which friction and gravity effects 
are taken into account in a simple way, through the parameter p. The value 
of the parameter p decides, whether the grain will stop on the site, or roll 
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further down the slope. 

The two threshold ricepile model introduced in |P, 0| is defined on a 
onedimensional lattice of size L, with a wall at the zero position and open 
boundary at the other end. At the open boundary, particles are free to flow 
out of the system. 

As in the experimental set up, the system is driven by adding particles 
to the position one, at the closed end. Every time unit one particle is added. 
Two thresholds are defined in the model: the critical threshold Zc, which is 
the local condition for the onset of avalanche, and the gravity threshold Zg 
{zc < Zg). If the local slope 

Zi = hi — /ij+i, (2) 

(where hi, i = 1,2,3,...,L is a height profile of pile) is less then Zc, it is 
too low and the friction stops the grain movement. The grain resides on 
the position i. If Zi > Zg, the local slope is too high and the grain moving 
downslope is not allowed to stop at the i-th position. But, in the case, that 
Zc < Zi < Zg, the grain moves from i to i + 1 with probability p. Through the 
parameter p, effective friction is introduced into the model. All supercritical 
slopes can topple, with probability p, but for local slopes, which are too big 
{zi > Zg), the gravity becomes decisive and the site topples with probability 
p = 1.0. 

The dynamics of the ricepile model is as follows: 

1. Each avalanche starts at z = 1. If zi > Zc, the site one is activated and 
topples a particle to the next nearest position (two), with the probability p. 
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If even z\ > Zg] p = 1.0. 

2. Every particle, sliding from the position z to z + 1 activates three columns, 
namely i ~1, i and i + 1. The position « — 1 is activated, because it possibly 
can become supercritical , when removing a grain from the z - th column. 
Columns i and i + 1 are activated, because they are destabilized by sliding or 
stopping particles, respectively. In the next time step, all supercritical active 
sites topple a particle to the i + 1-st position with probability p {p = 1.0, if 

Zi > Zg). 

3. Step two is repeated, untill there are no active sites in the system, that 
means, untill the avalanche is not over. 

Changing slightly the local toppling rules, different versions of the ricepile 
model are defined |^ : 

a) The number of particles toppled from the position i is constant and inde- 
pendent on the supercritical local slope Zi - the model is called limited. 

b) The number of toppled particles is a function of the supercritical local 
slope Zi, the model is called unlimited. 

c) If the particle (or more particles) topples from the site i and moves only 
to the next nearest position i + the model is defined as local. 

d) The model is called nonlocal, if n toppled particles moving from the i-th 
site are added subsquently to n nearest downslope positions (one particle per 
site) i + 1, i + 2,..., i + n. 

Thus four different ricepile models are recognized: 
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1. local limited model (LLIM) 

2. local unlimited model (LUNLIM) 

3. nonlocal limited model (NLIM) 

4. nonlocal unlimited model (NUNLIM) 

Universality classes for the ricepile model were studied by Amaral and 
Lauritsen [§] . Their results show, that local models (LLIM, LUNLIM) belong 
to the wide universality class called local linear interface universality class 
(LLI class) [|, 13] . The authors also found, that the nonlocal toppling rules 
lead to two new universality classes, with a different set of critical exponents. 
But none of the universality classes is the one of the real ricepile. 

3 Absorbing model 



Our absorbing model [Tl| exhibits simplified, one threshold version of the 
ricepile model . The gravity threshold is removed and all supercritical 
active sites are alloved to topple with probability p < 1. There , therefore, 
persists small, but nonzero probability, that also the extremely large local 
slopes are possible. Physically this seems to be quite plausible. It is not 
probable, that in real piles of granular material there exists a strict gravity 
threshold. There is rather a continuous transition to the local slopes , which 
are already so large that they, when activated, almost always topple. 

We investigated numerically all four versions of the absorbing model: 
LLIM, LUNLIM, NLIM, NUNLIM. Several quantities were measured for all 
of the models: 
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1. Material transport as a ratio of the number of outgoing to ingoing particles 

and its dependence on the parameter p. 

2. Average material transport {J{p)) as a function of p. 

3. Avalanche size distribution for different parameter values p. Avalanche 
sizes are measured in terms od dissipated potential energy, in accordance 
with the experiment 0. 

4. Changes in the pile profile, due to changes in the parameter p value. 

If the probability parameter p changes slowly in the interval (0, 1), the 
model typically passes through different dynamical regimes: 

i) isolating, in which all particles are absorbed in the system and none of 
them reaches the open boundary; 

ii) partialy conductive, in which the pile profile grows up as a bulk, because 
a certain fraction of the particles, depending on p, is absorbed in the system 
(absorbing properties); 

iii) and totaly conductive, when the number of ingoing and outgoing particles 
is balanced. 

4 Universality classes for the absorbing 
model 
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4.1 LLIM 

Dynamical properties of the local limited absorbing model are in details 
described in []TT| . Here I only briefly list the main results. 

Local limited toppling rules are defined as follows: active supercritical 
site topples one particle to the next nearest position with the probability p. 
Looking at the average material transport {J{p)), three dynamical regimes 
of the LLIM model are recognized {Fig. la): 

a) For < p < p ; p' ~ 0.53865, the system is completely isolating. The 
average transport {J{p)) is zero. For p close to zero, almost all ingoing par- 
ticles are absorbed. The avalanches die out soon , their size is exponentially 
bounded. 

Close to the first phase transition point p , the steepness of the pile is 
still high enough to say, that the local slopes are almost everywhere higher 
than the critical threshold Zc- This is the reason that the spreading of active 
sites in time is practically determined by the probability p; the same way as 
it is in the percolation process. In the space - time coordinate system, we 
have therefore a picture of directed percolation with three descendants and 
absorbing boundary |T3] . p' is thus simply the critical percolation threshold. 



Close to the percolation threshold p , the average transport J{p) scales with 
p as 

{J(p))-j'^(p-p/ (4) 
6' = 0.9 ±0.01 

where j' is the current flowing due to the finite size of the system. 



b) For the probability interval p < p < pc', Pc ~ 0.7185, the system is 
partialy conductive, with constant average slope. That means, the height 
profile grows as a bulk with velocity v{p). Fluctuations of transport J{p) 
exhibit long range correlations. 

Above the percolation threshold p , the percolation picture breaks down. 
The subcritical, absorbing states are randomly distributed throughout the 
system and the avalanche can stop anywhere. As p ^ Pc, the long range 
correlations in transport fluctuations are destroyed, and the region of small 
local slopes spans the whole system. The pile stops to grow and aX p = Pc 
it is pinned at the position i = L. Critical point pc is thus understood as 
the depinning transition point. Close to the depinning critical point, average 
transport scales as 

l-{J{p))^\p-Pc\' (5) 
6 = 0.9 ±0.01 

c) In the interval {pc, 1) the system is completely conductive. Transport 
fluctuations are of white noise type and the average transport J{p) = 1.0. In 
the dynamical regimes b) and c) the system is in the SOC state, having power 
law distribution of avalanche sizes (0) with critical exponent r = 1.57 ±0.05. 

4.2 LUNLIM 

Local unlimited toppling rules are in the absorbing model defined as follows: 
In order to get a realistic profile of the pile, each supercritical active site 
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topples k, k = mt(j^), grains to the next nearest downslope position with 
probabihty p. This way one gets smooth profile without cavities [Fig. 2a). 

Numerical investigations of {J{p)) reveals, that only two dynamical 
regimes are clearly recognized {Fig. lb): the pile is either completely isolating, 
{{J{p)) = 0), or completely conductive = 1.0). Partialy conductive 

dynamical regime is missing. 

a) For < p < pc, Pc ^ 0.6995 {Fig. lb), the system is completely isolating. 
From the definition of local toppling rules it is clear, that absorbing states 
{zi < Zc) are easily created even for very small values of the parameter p. It 
means, that the percolation picture in the space - time coordinate system is 
not correct in the case of the LUNLIM model. Near transition point Pc the 
average transport J{p) scales with p as {Fig.3) 

J{p)^{p-p,Y (6) 

6 = 1.93 ±0.07 

b) In the second dynamical regime {pc < p < 1.0) the system is completely 
conductive, with pile profile pinned at i = L {Fig. 2a). J{p) as a function 
of time exhibits white noise features {Fig. 4). Avalanche size distribution is 
critical (|I]) with critical power law exponent r = 1.54 ± 0.02 {Fig. 5a). 

4.3 NLIM 

Nonlocal limited toppling rule means, that the supercritical active site topple, 
with probability p, N particles to the nearest downslope positions. The 
nonlocal limited toppling rules preserve three dynamical regimes, the same 
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way as it is in the LLIM case. In Fig.lc, isolating, partialy conductive and 
totaly conductive regimes are recognized. 

a) For < p < p , p' ~ 0.267 (Fig.lc), the pile is in the isolating regime. To 
understand the nature of the first phase transition point p , the percolation 
picture in the time - space coordinate system is still usefull. But now, the 
number of descendants is, in principle, greater than three. That is the reason 
for the fact, that the percolation threshold is shifted to the lower parameter 
values as one can see when comparing Fig.la and Fig.lc, e. g. (p'nlim < 
p'llim)- III the model studied here, the number of particles toppling from the 
activated supercritical site is 4. Five sites are thus activated by every toppling 
from the position i, namely z — 1, z, z + 1, z + 2, z + 3. There are therefore 
five descendant sites in the directed percolation with absorbing boundary 



11] . In order to estimate p with greater accuracy, systematic studies of the 
dependence of percolation threshold on the number of descendant sites are 
necessary. 

The average transport near the percolation threshold p scales as (Fig.Qa): 

J{p)^{p-pY' (7) 

with the critical exponent 

5' = 1.18 ±0.04 

b) In the interval p < p < pc, the pile grows up with constant velocity v{p), 
maintaining the global slope on a constant value for a constant probability 
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parameter p. Transport J{p) as a function of time shows long range cor- 
relations, on the contrary to the totaly conductive regime, where it has a 
character of white noise {Fig. 7). 

c) Depinning transition occurs at Pc ~ 0.365 {Fig. 6b). Average transport 
scales with p close to the critical point as : 

l-{J{p))^\p-Pc\' (8) 
S = 1.12 ±0.06 

For the probability interval pc < p < 1, the profile of pile is pinned at 
i = L {Fig.2b), and the average transport {J{p)) = 1.0 {Fig.lc). Fig.5b 
demonstrates the avalanche size distribution in a case of partialy conductive 
and conductive dynamical regimes. In both cases the dynamics of the pile 
is self organized critical, which is demonstrated by the critical , power law 
scaling (|l|). The critical exponent r = 1.35 ± 0.05. 

4.4 NUNLIM 

The nonlocal unlimited toppling rule is defined as follows: N{zi) particles are 
released (with probability p) from the activated supercritical position and are 
added to N{zi) nearest downslope positions. 

The nonlocal unlimited version of the absorbing model shows completely 
different behaviour. First, no distinct dynamical regimes are recognized. The 
pile is completely conductive already for p close to zero as can be seen from 
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Fig. 8a. Pile profile is pinned at z = L {Fig. 2c). Avalanche size distribution 
shows the critical scaling (|l|) with critical exponent r = 1.51 ±0.05 {Fig.be). 

5 Discussion and conclusion 

Probability density function (|l|) scales with the system size as 

L) = L-^9{^) (9) 

with (3 = Dr. For the LLIM, LUNLIM and NUNLIM absorbing models the 
best data collapse has been found for D = 2.24, what indicats, that these 
models belong to the same universality class, called LLI universality class 

On the contrary, for the NLIM absorbing model the best data collapse 
was found for D = 1.55. The critical exponents r and D are different from 
that of LLI class and defines a new universality class to which belongs also 
the NLIM version of the two threshold ricepile model [§] . 

The reason of lowering the r exponent of the NLIM model in comparison 
with the LLIM model is as follows: The average slope of the NLIM and the 
LLIM pile is simillar. For example for p = 0.8, the average slope of the NLIM 
pile is 67.87° and for the LLIM pile 60.53°. Due to the nonlocal toppling rules 
in the NLIM model, more columns are perturbed and thus the probability of 
greater avalanches is enhanced. Therefore the exponent r is lowered. 

The same argument could be used in the case of the NUNLIM and the 
LUNLIM model. But here the situation is different. The average slope 
changes significantly with changes in toppling rules from local to nonlocal. 
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For example if p = 0.8, the average slope of the LUNLIM model is 82.22" 
and that of the NUNLIM model equals 53.6°. Because the number of toppled 
particles is proportional to the slope in the unlimited model, it seems, that the 
relatively small average slope of NUNLIM pile leads to relatively few particles 
released on average in one toppling. This fact should enhance the probability 
of small avalanches and the avalanche size distribution function should have 
T exponent greater than 1.55. This really happens for the NUNLIM two 
threshold ricepile model |^ . In this model the propability of big local slope 
decreases exponentially with Zi. But looking at {Fig. 2c) one can see, that it 
is not an exception to have a big local slopes in the NUNLIM absorbing pile. 
During a toppling event, the site with big local slope releases a number of 
particles (proportional to the local slope), which disturb a lot of downslope 
columns. This efect increases the probability of large avalanches. It seems, 
that in the case of the absorbing model the two described effects balance 
each other and thus the exponent r remains untouched by changed toppling 
rules. This is different from the NUNLIM two threshold model P]. Here 
the first effect is decisive and the model belongs to a new universality class 
(r = 1.63). 

Another question, which should be discussed, is the nonexistence of par- 
tialy conductive regime of the LUNLIM model. First, the model is local. 
That means, in every toppling, only three columns are activated by each 
toppled particle. Therefore, the probability of avalanches, having a chance 
to reache the end of the system and thus to transport a material, doesn't 
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increase due to more activated sites by every toppled particle. As it has been 
already told, in the LLIM model, there are no absorbing {zi < Zc) states in 
the system for the isolating dynamical regime. This is not the case of the 
LUNLIM model. Absorbing states, therefore, exhibits another obstacle for 
bigger avalanches to develop and transport the particles. The existence of 
absorbing states, even for a small parameter values, destroys the percolation 
picture of the spreading of active sites and this is also the reason for nonex- 
isting critical percolation probability p and only the transition to completely 
conductive dynamical regime is present. 

Last, some words should be told about the pile slopes in all of the three 
dynamical regimes. In the isolating regime average transport {J{p)) — 0. 
All the added particles are absorbed in the system. Moreover, the avalanche 
sizes in this regime are exponentialy bounded. That means, majority of the 
particles is absorbed on the first few columns of the pile. Therefore, if the 
driving time t tends to infinity, the average slope of the pile grows to infinity. 
That in consequence means, that also the local slopes become arbitrarily 
large. 

In the partialy conductive regime, constant amount of particles, depend- 
ing on the parameter p, is absorbed in the system. As t ^ oo, average slope 
of the pile remains constant; pile grows up as a bulk. The local slopes are 
finite, except of z{L) (see (2)), which tends to infinity. 

In the conductive regime {J{p)) — 1.0. The average slope of the pile is 
constant depending only on the parameter p. All local slopes are finite in 
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this regime. 

In conclusion, we have studied numerically the LLIM, LUNLIM, NLIM 
and the NUNLIM absorbing models. We have found, that the models belong 
to two different universality classes, characterized by different critical expo- 
nents. Both of the universality classes are different from the one of real pile 
of rice. We have studied the transport properties of all of the models and 
found phase transitions between different dynamical regimes. We state, that 
the dynamics of LLIM and NLIM model is directly mapped to the dynamics 
of directed percolation process at the absorbing boundary JTSf for the defined 
interval of parameter p. 

This work has been supported by the grant of VEGA number 2/6018/99. 
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6 Figure captions 

Fig.l 

Average transport ratio of particles through the system as a funcion of 
the parameter p. 

(a) LLIM model: Three different dynamical regimes are recognized: isolating, 
p e (0,0.53865); partialy conductive, p e (0.53865,0.7185) and conductive, 
p e (0.7185,1.0). 

(6) LUNLIM model: Only two different dynamical regimes are recognized: 
isolating, p G (0,0.6995) and conductive, p e (0.6995, 1.0). 
(c) NLIM model: Again, three different dynamical regimes are depicted: 
isolating, p e (0,0.267); partialy conductive, p e (0.267,0.365) and totaly 
conductive, p e (0.365, 1.0). 

Fig.2 

Pile profiles of the LUNLIM {Fig. 2a), NLIM {Fig. 2h) and NUNLIM 
{Fig.2c) absorbing models for different values of the probability parame- 
ter p. Notice, that in the totaly conductive regime, the pile profile is pinned 
at i = L (o), (c). The pile is growing as a bulk with velocity v{p) in the 
case of the partialy conductive regime and is pinned in the totaly conductive 
regime (6). 

Fig.3 

Ln-ln plot of the average transport as a function of the distance from the 
critical point Pc in the LUNLIM absorbing model, e— {p — pc)- We find, that 
the best scaling is obtained for p^ = 0.6995. 
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Fig.4 

Transport J{p) as a function of time (in iterations) in the conductive 
regime of the LUNLIM model for two different values of p. For p > pc, 
Pc = 0.6995, white noise is observed. For p < pc all particles are absorbed 
and therefore there are no fluctuations. 

Fig.5 

Ln -In plots of the power law parts of the avalanche size distributions 
(unnormalized). The critical exponent r of the LLIM LUNLIM (a) 
and NUNLIM (c) models is r = 1.55. The NLIM model (b) belongs to the 
different universality class with r = 1.35. The system size in Figs.5b,c is 
L = 300 and in Fig.5a L = 500. In the NLIM model (b) the bump shift 
with system size has been numericaly tested. With growing system size, the 
bump shifts to the higher values of s, indicating thus SOC. 

Fig.6 

(a) Ln-ln plot of the average transport as a function of the distance from the 
critical point p in the NLIM absorbing model, e = (p — Pc)- We find, that 
the best scaling is obtained for pc = 0.276. 

(b) Ln-ln plot of the average transport as a function of the distance from the 
critical point pc in the NLIM absorbing model, e = {pc~ p)- We find, that 
the best scahng is obtained for p^ = 0.3441. 

Fig.7 

Transport J{p) as a function of time (in iterations) in the partialy con- 
ductive regime {p = 0.339) and totaly conductive regime {p = 0.4) of the 
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NLIM absorbing model. For p > pc, Pc = 0.365, white noise is observed. 
For p < p < Pc the character of fluctuations is different. Long range corre- 
lations (reminiscent of Brownian motion) are observable in the time signal. 
For p < p all particles are absorbed and therefore there are no fluctuations. 

Fig.8 

Transport J{p) as a function of time (in iterations) in the NUNLIM model 
for three different values of the parameter p. The pile is totaly conductive in 
wide range of the parameter p and the signal has a character of white noise, 
with fluctuations depending on p. 
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